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Abstract 

We formulate low-frequency charge-fluctuation forces between charged cylinders 
parallel or skewed in salt solution: forces from dipolar van der Waals fluctuations 
and those from the correlated monopolar fluctuations of mobile ions. At high salt 
concentrations forces are exponentially screened. In low-salt solutions dipolar en- 
ergies go as -R -5 or i?~ 4 ; monopolar energies vary as i? -1 or lni?, where R is the 
minimal separation between cylinders. However, pairwise summability of rod-rod 
forces is easily violated in low-salt conditions. Perhaps the most important result 
is not the derivation of pair potentials but rather the demonstration that some of 
these expressions may not be used for the very problems that originally motivated 
their derivation. 
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Models of assembly by rod-like particles often require closed-form expressions for 
molecular interaction. In addition to brute steric forces, hydration forces, and electrostatic 
double layer interactions H between charged rods, there are interactions of correlated 
charge fluctuations 0. These fluctuations can be either the dipolar events of traditional 
van der Waals forces or monopolar charge fluctuations from transient changes in the 
density of mobile ions around a charged particle |§. 

Between like particles, these correlated-fluctuation forces are attractive. There is 
always a question, though, whether one is allowed to compute these attractive forces in 
molecular arrays as though they were the sum of individual rod-rod interactions. In this 
paper we derive the interaction between pairs of cylindrical rods - parallel or skewed - in a 
salt solution. We find that the conditions for the validity of pairwise additivity can be quite 
restrictive. Without the screening effect of the intervening salt the counterion correlation 
forces come from collective ionic fluctuations in the array that can not be decomposed 
into separate pairwise additive contributions. The range of validity of different theories 
||, that formalize counterion correlation forces is practically limited. 

We use the Lifshitz-Pitaevskii approach |3], [5| that begins with an artificial but 
easily formulated interaction between arrays of cylindrical rod-like molecules. This ap- 
proach implicitly reveals any dilute limit in which the pair potential between rods can be 
extracted from the array-array interaction. Specifically, we examine the van der Waals 
attraction between two like anisotropic media, L and R, containing parallel cylindrical 
particles embedded in a medium m (Figure 1). L and R are separated by an isotropic 
region of salt solution devoid of cylinders. 

The anisotropic regions are composites of parallel cylindrical rods of radius a at 
volume fraction v = N na 2 ; N is the cross-sectional density of N rods per unit area. The 
cylinders c have anisotropic intrinsic susceptibilities e c ± and e|, perpendicular and parallel 
to the rod axis. For each region r =L, m, R the effective dielectric and ionic properties of 
the composites can be written in terms of their local values. In this construction the axis 
parallel to the rod in region R is rotated about axis z to create an angle 9 with respect 
to the rod x-axis in L. 

Because ionic fluctuations are slow [Q, we consider only 'zero-frequency' van der 
Waals interactions, describing thermodynamic as opposed to quantum fluctuations. The 
wave equations in such media follow 



where k is the same as the Debye screening constant except for the missing dielectric 
constant e in the denominator. The dielectric susceptibility tensor e is taken only in the 
limit of zero frequency. The susceptibility tensors e r for the composite regions r = L and 
R have diagonal elements e r x , e T y , e r z , and with the axes x and y of R rotated by an angle 
9. For small enough volume fractions v, the components of the unrotated tensors for the 
composite media are H e L and e R 
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A II = J1— - A 1 = ^-^ (4) 



where e m is the isotropic susceptibility in region m and between cylinders in L and R; 
superscript c refers to the response of the cylinder material. 

The quantities k 2 in each region r = L,m, R depend on volume averages of ionic 
strengths n = Z 2 nz where nz is the average number- density of ions of valence Z. 
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For region m, k^ = 47r "" e , the Debye constant being defined in a standard way as 
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For regions L and R, we distinguish the densities of mobile ions n c and n m in 
the cylinders and in the intervening space respectively. The mobile ions associated with 
the cylinders model the tightly associated counterions where n c pertains to the average 
density of mobile charges within the cylinders. The volume average ni and n R over the 
whole composite is n L = tir = n c v + n m (l — v) = n m + v{n c — n m ), which leads to 

kl = k 2 R = (n m + v(n c - n m )) = k 2 m (l + v ^- ^ . (5) 

Because solutions of the wave equation require continuity of eV0 perpendicular to the 
interfaces at z = and I, it is the perpendicular response of e L and that couples with 
k\ = k\ to create the ionic screening lengths in both semi-infinite regions. For small 
enough volume fractions v 
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From this point on we will simply write k 2 (no subscript) for n Tl 
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The 'zero-frequency' van der Waals interaction free energy per unit area between 
the two anisotropic semi-infinite regions L and R across a slab m of thickness I is || 
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(1,0) = -y J j^y 2 ln(2?(g,Z)) = I # I HV{Q,l))QdQ (7) 



Here, the Q are radial wave vectors within the x, y plane; -0 is for angular integration 
over all directions in Q to take account of anisotropy. The secular determinant for the 
relevant modes can be written as 



V(Q, = 1- A Lm ^)A Rm (6 - V) x e - 2 ^^ 1 . (8) 
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The functions A£ OT (^>) and A^ m (6* — xji) necessarily go to zero as N,v — > 0; both have 
the form (i = L, R) 



SrnVW+K* ~ ejQ 2 (l + ( °-^) COS 2 if) + & 

A iM = V , V V _ X ' ' ± 0) 



For small enough values of v, the A im ((f)) = A im (<j), v) functions can usually be expanded 
to terms linear in cylinder volume fraction v. To lowest order in density, the interaction 
energy goes as iV 2 . In this dilute limit, rods in the two media interact pairwise across the 
gap I. 

Extraction of the pairwise interaction potential follows one of two procedures, 
depending whether or not the rods in L and R are parallel (8 = 0) or are skewed, to 
connect the per-unit-area interaction Q{1, 8) between planar regions of embedded cylinders 
with the pair interaction potential g(l, 8) between skewed cylinders, or the pair interaction 
energy per unit length g(l,8 = 0) between parallel cylinders. A sufficient condition for 
this connection is that Q(l,8) be accurately expressed by the first, quadratic term in a 
series expansion in density iV or volume fraction v. 

The Pitaevskii ansatz || can be applied in two forms. For skewed cylinders at a 
mutual angle 8, the connection is 



= N 2 sin 8 g(R, 8) + 0(N 3 ) (10) 

l=R 
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while for parallel cylinders the connection is an Abel transform 

d 2 g(i,d = o) 



lim 
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N 2 J g(^/R 2 + y 2 , 8 = 0) dy + 0(N 3 ). (11) 



n^o dl 2 

In both cases R is the minimal separation between the cylinders. 

Expanding Q(l, 8) to the second order in v one obtains 
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For skewed cylinders this leads to 
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for parallel cylinders through the inverse Abel transform one remains with 

kT r°° 

g (R, 9 = 0) = -2— K 5 (7ia 2 ) 2 / p 4 dp e- 2KR ^K (2 K Rp)f(p, 9 = 0) (15) 
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K (x) is the modified Bessel cylindrical function of order 0. 



Integrating over the azimuthal directions -0 and taking the second derivative 
d g Jp^ , reveals three classes of interaction depending on the type of charge fluctuation: 
dipole-dipole, dipole-monopole, and monopole-monopole: 
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We examine the results for g(R, 9) and g(R, 9 = 0) in several limiting cases and 
then ask when one is permitted to use expressions that pre-suppose pairwise additivity of 
rod-rod interactions. 



In the limit of vanishing salt concentration in the medium, 2kR — > 0, 
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When all salt concentrations equal zero, interactions come only from differences 
in dielectric susceptibility. All monopolar terms vanish to leave only dipolar terms 0. 

In the opposite case of strong ionic screening, 2kR 3> 1, the only terms surviving 
are those of longest range that nevertheless decay exponentially. 
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When all e's are equal but there are differences in ionic densities n c and n m , only monopo- 
lar terms endure, 
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There is a curious further limit - unscreened monopole-monopole correlated fluc- 
tuations - where n m — > but cylinders still carry counterions n c , 
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The n m = form for parallel rods has been obtained analytically [|J before and 
checked against computer simulation [[|. Our result Eq.E3 however indicates that only 
the linearized form of the Oosawa potential is pairwise additive in an array, i.e. there is 
no leveling off of the correlation attraction as the concentration of mobile counterions, n c 
in our model, increases. 



That Eqs. g3j and g4| indeed correspond to monopolar fluctuations can be es- 
tablished by recalling |7] monopolar charge fluctuation forces between small spheres of 

radius a at separation R that go as ~ kT (j|J e~ 2nR . Assuming that this ion correlation 

potential sums pairwise between different parts of the rods in either skewed or parallel 
configuration, one obtains the forms of Eqs. B3| and p4. 



When can one use pairwise rod-rod interactions in solutions and in arrays? 



The central requirement of the Pitaevskii construction is to create composite 
media whose dielectric properties are linear in the number density of interacting particles. 
The interaction between composites then goes as the square of the number densities. With 
salt solutions, requisite dependence on the square of densities requires more than linear 
dielectric response. 



Expanding the dispersion relation T>(Q, I) (Eq. in powers of cross sectional 
density of the cylindrical array, we see that the expansion 
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are strictly satisfied. 



Whenever salt concentration n m — > while the concentration of the mobile 
charges within the cylinders n c remains finite, expansion in N must always violate the 
first inequality. Bluntly, limiting forms such as Eqs. ^ and cannot be used in problems 
where there are more than two molecules. Their very range carries within it the source 
of their self-invalidation as pairwise-additive interactions. 



Traditional van der Waals forces, due to differences in dielectric susceptibility, 
will often not be pairwise additive. For low-frequency fluctuations, the second inequality 
in Eq. |26| can easily be violated 0. Finite values of the medium salt concentration 
always create an exponentially damped contribution when 2n m R ^> 1 . Though screening 
severely weakens forces, it does act to make pairwise summation more accurate. 



Charge fluctuation forces in condensed systems can be qualitatively non-additive. 
This feature imposes questions that are sometimes neglected in analyses of colloidal and 



6 



macromolecular systems. When addressing correlated ionic fluctuations it appears to be 
necessary to think in terms of collective excitations that extend over all molecules in an 
array. 
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Figure 1: Array-interaction geometry. Region L is confined to z < 0, region R to z > I; 
isotropic region m is the slab in between. The cylindrical array in R is rotated about the 
z axis by an angle 9 with respect to L. Not only the dielectric constants of the composite 
media R and L but also the dielectric constants of the cylindrical rods themselves are 
anisotropic, ejj, e c ± , n c refer to the cylinder material. 
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